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Abstract 

We study the spreading of an infection within an SIS epidemiolog- 
ical model on a network. Susceptible agents are given the opportunity 
of breaking their links with infected agents, and reconnecting those 
links with the rest of the population. Thus, the network coevolves 
with the population as the infection progresses. We show that a mod- 
erate reconnection frequency is enough to completely suppress the 
infection. A partial, rather weak isolation of infected agents suffices 
to eliminate the endemic state. 



1 Introduction 

Among the many potential applications of agent-based models whose pattern 
of interactions is represented by a network, a broad class is constituted by 
systems where that pattern is not a static structure, but evolves in response to 
the changing state of the agents. Generally, the evolution of the interaction 
network and the dynamics of individual agents occur over different time 
scales. For instance, in learning processes -such a those implemented in 
real and artificial neural networks |1] connections change adaptively over 
scales that are large as compared with the internal dynamics of agents. At 
the opposite limit, in models of network growth, the pattern evolves in the 
absence of any dynamics related to agents at the network nodes [2]. 

When, on the other hand, the dynamical time scales of a population 
of agents and its interaction network are comparable, we can speak about 
their coevolution. Consider, for instance, agents involved in a cooperation- 
defection game, such as in the prisoner's dilemma [3J, where defected agents 
are given the opportunity of breaking the connection with their defectors, in 
such a way that no further interaction is possible between them. Or, in a 
process of opinion diffusion, that agents which do not succeed at reaching an 

1 This is a written version of a talk to be given at the International School on Com- 
plexity: Course on Statistical Physics of Social Dynamics: Opinions, Semiotic Dynamics, 
and Language (Erice, 13-20 July, 2007). 
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agreement may elude mutual contacts in the future [HElElEj. In such cases, 
a realistic assumption would be that the breaking of interaction links occurs 
just after a few, or even one of, such events. 

In this talk, we explore a model of agent-network coevolution in a popula- 
tion where an infection is spreading, such that non-infected agents are given 
the possibility of avoiding contact with their infected partners -perhaps in 
response to risk perception [8]. The model is based on an SIS epidemiological 
process, where each agent can be susceptible (S) or infected (I). We recall 
that in the standard SIS process each infected agent recovers and becomes 
susceptible at a fixed rate, say, with probability 7 per unit time. Susceptible 
agents, in turn, become infected by contagion from infected agents, at a rate 
proportional to the infection probability per unit time, p, and to the fraction 
of infected agents, raj. Within a mean- field description of the standard SIS 
model, the fraction of infected agents obeys 



where ns = 1 — n\ is the fraction of susceptible agents. In this description, 
for asymptotically long times, nj vanishes if p < 7. Therefore, the infection 
is suppressed as time elapses. If, on the other hand, p > 7, the fraction of 
infected agents approaches a finite value n} = 1 — 7/p > 0, and the infection 
is endemic. The transition between these two regimes occurs through a 
transcritical bifurcation. 

We introduce in the following an implementation of the SIS model on a 
network, with agents occupying the nodes and contagion taking place along 
the links. A susceptible agent can become infected only if it is connected to 
an infected agent. Moreover, susceptible agents have the opportunity, with 
a certain probability, of breaking their connection with an infected partner 
before contagion takes place, and reconnect the broken link with any other 
agent in the population. We show that these reconnection events are able to 
control infection spreading, even to the point of completely suppressing the 
infection as the reconnection probability grows. This suppression is achieved 
with a moderate decrease in the number of network connections per infected 
agent. 

2 Model and mean-field formulation 

The present model is a variation of the model introduced by Gross et al. [9] . 
Consider a population of iV agents at the nodes of a network with M links. 
The average number of neighbours per agent is z = 2M/N. At a given time, 



hi = — 7711 + pniiis, 
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each agent can be in one of two states: susceptible (S-agent) or infected (I- 
agent). Initially, the M links are distributed at random over the population, 
and there is a certain fraction of agents in each state. 

Each evolution step is divided into two sub-steps. In the first sub-step, 
an agent is chosen at random from the whole population. If it is infected, it 
recovers with probability 7, and becomes susceptible. In the second sub-step, 
a pair of linked agents is chosen at random. If both agents are susceptible 
or infected, nothing happens. Otherwise, the S-agent is given the chance 
to break the link with its infected neighbour and reconnect it with another 
agent, taken at random from the remaining of the population. This rewiring 
happens with probability r. Finally, if the rewiring has not occurred, the 
S-agent becomes infected with probability A. 

A time unit contains N evolution steps. Thus, each I-agent has a prob- 
ability 7 per time unit of becoming susceptible, so that the mean duration 
of the infection period is 7 -1 . Moreover, each S-agent in contact with an 
I-agent becomes in turn infected with probability 2(1 — r)A per time unit. 

2.1 Mean- field equations 

Although, clearly, the model is defined bearing in mind its implementation as 
a numerical simulation, a formulation in terms of differential equations, from 
mean-field arguments, turns out to give an essentially correct description of 
the system. To define our mean-field formulation, we first introduce a suitable 
set of variables. We call Ni and iVg the number of infected and susceptible 
agents, respectively. They satisfy A^i + A^s = N. Moreover, we denote by Mn, 
Mis, and Mss, respectively, the number of network links joining two infected 
agents (II), an infected agent and a susceptible agent (IS), and two susceptible 
agents (SS). Since the total number of links is preserved by the dynamical 
rules, we have Mn + Mis + Mss = M constant at all times. Differential 
equations will be formulated for the fractions n\ = Nj/N, mn = Mn/M, and 
mis — Mis/M. The conservation of the number of agents and links implies 
ns = N s /N = 1 — n\ and m S s = M S s/M = 1 — mn — mis- 

The evolution of the above defined quantities is given by the events of 
infection and recovery, and by the reconnection of links. When, for instance, 
an I-agent becomes susceptible, there is not only a decay in the fraction of 

I- agents, but also a change in the fractions mn and mis. In fact, the links 
joining the recovered agent with I and S-agents pass, respectively, from the 

II- type to the IS-type, and from the IS-type to the SS-type. A symmetric 
situation occurs when an S-agent becomes infected. The number of links 
of each type associated to a given agent is calculated using mean-field-like 
averages. For instance, the number of Il-links associated to an I-agent is 
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estimated as 2Mu/N\ = zmn/ni, where z = 2M/N is the overall average 
number of links per agent. Similarly, the number of IS-links associated to 
an I-agent is taken to be Mis/Ni = zmis/2rii. Reconnection events, in turn, 
can change an IS-link into an SS-link. The probability for this change is 
proportional to the fraction of S-agents in the whole population and to the 
reconnection probability r. Finally, it has to be taken into account that both 
infection and reconnection only occur when an IS-link is selected, so that 
their probabilities are proportional to mis- 

Putting all these considerations together, the evolution equations for m, 
mn, and mis read 

n[ = —m + Amis, 

m' u = -2m n + Am| s / (1 - m), ,^ 
m' ls = 2m n - mis - f(l - ™i)™is 

+Ami S (2-2m n -3mi S )/(l -n{), 

with A = (1 — t)A/7 and r = 2r/z , y. Primes indicate differentiation with 
respect to the rescaled time t' = yt. In this formulation, thus, the inverse 
recovery probability 7 -1 fixes an overall time scale, and the infection and 
reconnection probabilities are accordingly redefined as A and f, which are 
the only two parameters left. Note that A depends on both A and r, and 
that r incorporates the only network-specific parameter, namely, the average 
connectivity per agent z = 2M/N . 

In the absence of reconnection events, r = 0, and provided that we put 
Tnii = Hi an d m is — 2njns = 2m(l — rii), the three lines in Eq. (TSJ) collapse 
into 

n[ = —rii + 2Ani(l — rii). (3) 

Rewriting this equation in terms of the non- normalized parameters, we reob- 
tain the mean-field equation (TjQ) for the SIS model, with p = 2A. The factor 
of 2 relating the infection frequencies p and A originates in the fact that, in 
our implementation of the SIS model on a network, each agent is (on the 
average) chosen twice per unit time to possibly become infected. 

2.2 Infection level at equilibrium 

We focus the attention on the equilibrium solutions of Eqs. (j2J), which are the 
candidates to represent the infection level and network structure at asymp- 
totically long times. First, we consider the stationary values of the fraction 
of infected agents. The analysis is restricted to the case of A > 1/2 which, in 
the absence of reconnection events, corresponds an endemic infection where 
a finite fraction of the population is infected at all times, ni > for t —>■ 00. 
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At the fixed points of the dynamical equations ([2]), the equilibrium frac- 
tion of links mjj and rajfg are related to the equilibrium fraction of infected 
agents n\ as 

m * = -_^L m * = 5. ( 4 ) 

11 2A(1 - n\) A V ; 

In turn, nj 1 satisfies 

= nj 2A - 1 - f + (3f - 2A)nJ - 3fnf + fnf . (5) 

This polynomial equation has four solutions. One of them tends to infinity 
for f — > 0, and remains real and larger than one for any positive f. Since 
meaningful solutions to our problem must verify n{ < 1, we disregard this 
solution from now on. 

The trivial solution n[ ^ = exists for any value of the normalized in- 
fectivity A and of the normalized reconnection probability f. For a given A, 
its stability depends on f. As discussed in more detail below, = is 
unstable for small f and becomes stable as f grows. The other two solutions 
read 



,(1.2) 

with 




n{ ' ' = l — \ j — cos — =F V 3 sin — (6) 



'32A 3 

a = arctan \l 1 (7) 

27f 

(0 < a < 7r/2). These two solutions are real for 32A 3 > 27f. Otherwise, 
they are complex conjugate numbers. The solution n[ , with the minus 
sign in the right-hand side of Eq. (jSJ), approaches 1 — (2A) _1 for f — *■ 0. 
Thus, it represents the expected fraction of infected agents in the absence of 
reconnection. When it is real, it satisfies rij < 1, and it is stable as long 
as it remains positive. Consequently, along with the trivial solution, 

(2) 

is another meaningful equilibrium solution to our problem. Finally, n\ is 
negative and stable for small f. Depending on A, it can become positive as 
f grows but, at the same time, it becomes unstable. Therefore, it does not 
represent a meaningful solution. 

Figure [T] summarizes, in a bifurcation diagram, the behaviour of nf \ 

(2) 

and n\ as functions of the normalized reconnection probability f, for three 
representative values of the normalized infectivity A. In the three cases, we 
have A > 1/2, so that -as discussed above- a non-trivial meaningful solution 
does exist. Full and dotted lines represent, respectively, stable and unstable 
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Figure 1: Bifurcation diagram for the equilibrium fraction of infected agents 
n\ as a function of the normalized reconnection probability f , for three nor- 
malized infectivities A. Although only positive values of n\ are meaningful, 
an interval in the negative domain is also shown for completeness. Full and 
dotted lines represent, respectively, stable and unstable branches. For clarity, 
the solution rij = is plotted in the vicinity of the transcritical bifurcation 
only. 
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branches. For small infectivity (A = 0.6), the stable solution nj crosses 
and becomes negative and unstable, while becomes stable. This 
transcritical bifurcation takes place at f = 2 A — 1. As f grows further, 

(2) 

and the negative stable solution n\ approach each other, and collide when 
rij = = 1 — 3/4A. Beyond this tangent bifurcation, which takes place 
at r = 32A 3 /27, the two solutions are complex numbers. 

The situation is different for larger infection probabilities, as illustrated by 
Fig. [T]for A = 1. Now, for f = 0, np is large and, as f grows, it is the stable 
negative solution nj 2 ^ which first reaches n[ ^ . At the transcritical bifurcation 
at r = 2 A — 1, nf^ becomes positive and unstable, and nf^ becomes stable. 
The tangent bifurcation where rvp and nf^ collide and become complex, at 
f = 32A 3 /27, takes now place when these two solutions are positive. As a 
consequence, there is an interval of normalized reconnection probabilities, 
between the two bifurcations, where the system is bistable: both nf^ 1 and 
7ij are stable meaningful solutions to the problem. The asymptotic state 

(2) 

is selected by the initial condition for nj, and n\ stands at the boundary 
between the two attraction basins. 

The regimes of small and large infectivity are separated by the critical 
value A = 3/4 = 0.75, also shown in Fig. [TJ At this critical point, the 
transcritical and the tangent bifurcation collapse into a pitchfork bifurcation 
at f = 1/2. Here, the three equilibria collide simultaneously, and nj '* becomes 
stable, while the other two solutions become complex. 

A phase diagram of our system over the parameter plane (f, A) is shown 
in Fig. [21 The zones of endemic infection, where the fraction of infected 
agents at asymptotically long times is positive (ni — > n^), and of infection 
suppression (nj — > 0) are separated, for large A and f, by the bistability 
region, where both asymptotic behaviours can be obtained, depending on the 
initial condition. The three zones are limited by the lines of the transcritical 
bifurcation [A = (1 + r)/2, TC], where n| ^ = changes its stability, and of 
the tangent bifurcation [A = (27f/32) 1//3 , Tg] where nj 1 and nf 1 collide and 
become complex. These two lines are tangent to each other at the "triple 
point" (1/2,3/4), where the bistability region disappears, and the system 
undergoes a pitchfork bifurcation [P]. For smaller A and f bistability is no 
more possible, and the zones of infection persistence and suppression are 
separated by the transcritical line. The tangent bifurcation takes now place 
at negative values of nj 1 ^ and nf^ 1 (dotted line). 

The inset of Fig. [2] shows the same phase diagram in terms of the orig- 
inal, non- normalized reconnection and infection probabilities, r and A, for 
a recovery probability 7 = 0.01 over a network with an average of z — 10 
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Figure 2: Phase diagram in the (r, A)-plane, showing the regions of infection 
suppression {n\ — > 0) and persistence {n\ — > n^), and the intermediate 
bistability zone. Their boundaries are given by the transcritical (TC) and the 
tangent (Tg) bifurcation lines, which collapse into a pitchfork bifurcation (P) 
at (1/2,3/4). The dotted line is the continuation of the tangent bifurcation 
line in the zone where is negative. The inset shows the same phase 
diagram in terms of the non-normalized parameters r and A, for 7 = 0.01 
and z = 10. 
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neighbours per site. Note that the relation between non-normalized and nor- 
malized parameters is not a mere change of scale, because both r and A enter 
the definition of A. 

Let us summarize our results on the persistence or suppression of the 
infection in terms of the non- normalized parameters. First, for small in- 
fectivity, A < 7/2(1 — r), the infection is always suppressed. In this sit- 
uation, the infectivity is just too small to sustain a finite infected popu- 
lation. For larger infectivities, on the other hand, the infection can be- 
come established, depending on the reconnection probability r. In the range 
7/2(1 — r) < A < 37/4(1 — r), the infection is endemic if reconnections are 
infrequent, r < z(X — j/2)/(l + zX). Otherwise, for sufficiently frequent re- 
connections, the infection dies out. The transition between both situations 
is continuous in the fraction of infected agents, and occurs through a trans- 
critical bifurcation. For even larger infection probabilities, A > 37/4(1 — r), 
the regimes of persistence (low r) and suppression (large r) are separated by 
a bistability zone, where the infection persists or dies out depending on the 
initial fraction of infected agents. The bistability zone is limited by the trans- 
critical bifurcation quoted above and a tangent bifurcation at a reconnection 
probability given by the solution to 27j 2 r = 16(1 — r) 3 A 3 . The discontinuous 
nature of the tangent bifurcation implies that the endemic state present in the 
bistability zone disappears abruptly at the boundary, with a finite jump in 
the asymptotic fraction of infected agents, from ni = 1 — \/ z(l — r)X/3r > 
to zero. 

2.3 Number of neighbours of infected and susceptible 
agents 

The variables mn and mis characterize how the structure of the network is 
related to the state of the agents. Reconnection events favor the growth of 
the number of SS-links at the expense of IS-links. Thus, for r > 0, S-agents 
should asymptotically posses relatively large numbers of neighbours. The 
equilibrium values mjj and mjg as functions of the equilibrium fraction n{ of 
I-agents are given by Eqs. PJ. These equations show, as expected, that the 
fraction of links connecting I-agents with any other agent is proportional to 
the fraction of I-agents itself. 

In order to introduce quantities that define the connectivity of I-agents 
and S-agents independently of their respective fractions, we consider the av- 
erage number of neighbours per agent of each type. For I-agents, for instance, 
the average numbers of infected and susceptible neighbours are 2Mn/Ni and 
Mfs/Nj, respectively. The average connectivity of I-agents, Zi, is the sum of 
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Figure 3: Connectivity of infected and susceptible agents, z\ and zs, relative 
to the overall connectivity z = 2M/N, for three values of the normalized in- 
fectivity A, as functions of the normalized reconnection probability f. Only 
the values corresponding to meaningful stable solutions for the fraction of 
infected agents at plotted. The connectivity of infected agents is not plotted 
beyond the threshold of infection suppression. The vertical dashed line rep- 
resents the finite jump in z$ at the tangent bifurcation where the solution 
disappears. 
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these two quantities or, equivalently, 



zi 1 

- = —~ , 8 

z 2A(l-n*) 

which gives the ratio between z\ and the overall average connectivity per 
agent at equilibrium. With analogous arguments for S-agents, their average 
connectivity reads 

Z JL = _J ^ f9) 

z 1-n* 2A(1 -nf) 2 ' 

Due to the conservation of the total number of links, Z\ and zs are univocally 
related. This relation can be obtained from Eqs. (jSJ) and (jUJ) by eliminating 
n*j, which yields 

z s = z l [l + 2\(l-z 1 /z)]. (10) 

In order to describe the correlation between the structure and the state of 
the population it is however useful to analyze both z\ and z$ as functions of 
the relevant parameters. Figure [3] illustrates the behaviour of z\ and zs, as 
described in the following, for the infection probabilities A already considered 
in Fig. [TJ 

For n\ = n[ , which stands for the stable equilibrium solution for low 
reconnection probabilities, both z\ and zs approach z as f — > 0. As expected, 
in the absence of reconnection events, there is no difference in the number 
of neighbours of infected and susceptible agents. Also, as f grows from zero, 
we have z\ < z < z$- We verify that reconnection tends to increase the 
connectivity of S-agents at the expense of I-agents. 

The other solution relevant to the epidemiological process, n* = = 0, 
corresponds to a purely susceptible population. Accordingly, we find zs = z. 
Note also that Eq. (jSJ) predicts Z\ = z/2X, but this value is never realized 
due to the total absence of I-agents in this state. 

For A < 3/4, the fraction of I-agents decreases monotonically with r and 
vanishes continuously at the transcritical bifurcation -or, for A = 3/4, at 
the pitchfork bifurcation. The connectivity of S-agents is zs = z both at 
f = and at the bifurcation. For intermediate values of the reconnection 
probability z$ is larger than z and attains a maximum. This maximum, 
which at first sight may result to be surprising, can be easily explained. In 
fact, to sustain a value of zs larger than the overall average z, it is necessary 
to have I-agents with a relatively low number of neighbours. As the infection 
is progressively suppressed by reconnection, the number of I-agents decreases 
and, accordingly, their contribution to the average number of neighbours per 
agent becomes less significant. At the bifurcation and beyond, S-agents must 
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account for the whole average, so that z§ returns to its value for r = 0, i.e. 
z s = z. 

The connectivity of I-agents, in turn, is a monotonically decreasing func- 
tion of f, and reaches z\ = z/2X < z at the bifurcation. This implies that, 
even at the threshold of infection suppression, I-agents maintain a finite num- 
ber of neighbours within the population. 

For A > 3/4, again, the connectivity z$ associated with the solution 
initially increases with f, and attains a maximum. In the subsequent decay, 
however, it does not reach z$ = z. In fact, rij disappears through a tangent 
bifurcation when it is still positive, so that the jump in the infection level 
is discontinuous. At the bifurcation, we find z$ = 2z(2A + 3)/9 > z. The 
connectivity of I-agents decreases with f and, at the bifurcation, its value is 
independent of A: z\ — 2/3. 

From the viewpoint of the interplay of the epidemiological dynamics and 
the structure of the underlying network, the most interesting result of this 
section is the fact that the infection dies out even when infected agents keep 
a substantial connectivity with the rest of the population. In the cases il- 
lustrated in Fig. [H], for instance, infected agents preserve more than 60 % of 
their connections at the threshold where the infection level vanishes. In other 
words, reconnection needs not to completely isolate infected agents to sup- 
press the infection. A moderate, partial isolation of the infected population 
is enough to asymptotically inhibit the endemic state. 



3 Discussion and conclusion 

What mechanisms are at work when the infection is suppressed by reconnec- 
tion, even when the connectivity of infected agents remains fairly high? To 
advance an answer to this question, it helps to consider a simpler dynamical 
system for the fraction of I-agents: 

n[ = -m + [2A - f(l - ni) 2 H(l - rii). (11) 

The right-hand side of this equation is just a rearrangement of that of Eq. 
(JSJ). The equilibria of Eq. ( fill are thus identical to the equilibria for n\ 
in Eqs. (|2l) . Moreover, their stability properties are also the same as in our 
original system. It is important to understand, however, that fllip and §2§ are 
not equivalent: they merely share the same equilibrium behaviour in which 
regards the fraction of I-agents. 

We immediately see that Eq. (llip can be put in the form of the standard 
mean-field equation (JTJ for a SIS process if we introduce the effective infection 



12 



probability 



PcS = 7 [2A - f (1 - ni f} = 2(1 - r)A - — (1 - ni f. (12) 

In Eq. (JH) the threshold of infection suppression, where the trivial equilibrium 
changes its stability, is given by p = 7. Imposing this same condition to p c s, 
we find f = 2A — 1. But this is precisely the suppression threshold in the 
system with reconnections. Therefore, with respect to the stabilization of the 
trivial equilibrium, the system (j2J) is effectively equivalent to the standard SIS 
model with infectivity p e g. The transcritical bifurcation of Eqs. (J2J), where 
= becomes stable, can be interpreted as a kind of continuation for 
r 7^ of the transcritical bifurcation of the SIS model without reconnection 
events. 

The interpretation of the tangent bifurcation where the endemic state dis- 
appears at a positive value of rij , for A > 3/4, is less direct. It can however 
be argued that the presence of such a tangent bifurcation, together with the 
transcritical bifurcation which stabilizes the trivial equilibrium, constitute 
the generic critical behaviour expected for a SIS model like Eq. (|11|) . with 
an infection probability which depends on the density of I-agents: 

n[ = -ni + p(rii)ni(l -ni). (13) 

Besides the trivial equilibrium, this equation has fixed points at the solutions 
of 

p(nf) = (1 - nf)" 1 . (14) 

Figure H] illustrates graphically two representative situations. The dotted 
curve is the graph of the right-hand side of Eq. (1141) as a function of n\. If 
the graph of p{n\) has a single intersection with the dotted curve (A) and 
if, upon variation of parameters in the infection probability, the graph varies 
as indicated by the arrow, the intersection crosses n{ = and a transcritical 
bifurcation takes place. The standard SIS model, for which p is constant, is 
a special case within this situation. More generally, the graph of p{n\) may 
have two (B) or more intersections with the dotted curve. When parameters 
change, it is still possible than one of the intersections becomes involved 
in a transcritical bifurcation crossing n{ = 0, exactly as in situation A. 
Now, however, it may well be the case that two intersections approach each 
other, and eventually collapse and disappear, as in B. In this case, Eq. ( TT3|) 
undergoes a tangent bifurcation, as found to happen in our system 

In summary, the above discussion shows that the suppression of the en- 
demic state as a result of reconnection events can be reasonably understood 
in terms of the critical behaviour of a standard SIS model with an effective 
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Figure 4: Graphical solution of Eq. ( TT4"I) . The dotted curve represents the 
right-hand side of the equation, and full curves are possible graphs of the 
left-hand side. Dots stand at their intersections. The arrows illustrate how 
the graphs may change upon the variation of parameters, in the cases of a 
transcritical bifurcation (A) and of a tangent bifurcation (B). 
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infectivity, which depends on both the reconnection probability and on the 
fraction of infected agents. The transcritical bifurcation which stabilizes the 
state where the infection is completely inhibited is interpreted as a contin- 
uation of a similar transition in the absence of reconnection. In turn, the 
tangent bifurcation -which, for large infectivities, suppresses the infection as 
the reconnection probability grows- is a generic phenomenon in SIS models 
with density-dependent infectivity [TO] . 

*## 

We have studied a model for an epidemiological process in a population 
of agents on a network, where contagion can occur along the network links. 
The network coevolves with the population as the infection progresses: a 
susceptible agent can decide to break a link with an infected partner, and 
reconnect it with another agent. Our main result is that this reconnection 
mechanism, implemented with moderate frequency, can completely suppress 
the endemic state where an infected portion of the population persists at 
arbitrary long times. Suppression of the endemic state does not require full 
isolation of the infected population. On the contrary, it can be achieved while 
each infected agent preserves a substantial part of the links with the rest of 
the population. 

Whether these results are relevant to real epidemiological processes is a 
question beyond the present study. Our analysis should be regarded as an 
illustration of the kind of collective effects that may emerge from the coevo- 
lution of populations of dynamical elements and their interaction network. 
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